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Self-gravitating systems such as elliptical galaxies appear to have a constant specific entropy 
and obey a scaling law relating their potential energy to their mass. These properties can 
be interpreted as due to the physical processes involved in the formation of these structures. 
Dark matter haloes obtained through numerical simulations have also been found to obey a 
scaling law relating their potential energy to their mass with the same slope as ellipticals. 
Since the X-ray gas in clusters is weakly dissipative, we have checked the hypothesis that it 
might verify similar properties. 

We have analyzed ROSAT-PSPC images of 24 clusters, and also found that: 1) the Sersic law 
parameters (intensity, shape and scale) describing the X-ray gas emission are correlated two 
by two; 2) the hot gas in all these clusters roughly has the same specific entropy; 3) a scaling 
law linking the cluster potential energy to the mass of the X-ray gas is observed, with the 
same slope as for elliptical galaxies and dark matter haloes. 



1 Introduction 



The optical surface brightness profiles of elliptical galaxies can be fit by a Sersic law (Sersic 
1968): 

(1) 

characterized by three parameters: So (intensity), a (scaling) and v (shape). We have shown for 
a sample of 132 ellipticals belonging to three clusters that the Sersic parameters were correlated 
two by two, and that in the three-dimensional space defined by these three parameters they were 
located on a thin line (Fig. |l]). These properties have been interpreted as due to the fact that, 
to a first approximation, all these elliptical galaxies have the same specific entropy (entropy per 
unit mass) (Gerbal et al. 1997, Lima Neto et al. 1999, Marquez et al. 2000), and that a scaling 
law exists between the potential energy U and the mass M for these galaxies (Marquez et al. 
2001): U oc jVf 1 ' 72=t0 ' 03 . Each of these relations defines a two-manifold in the [ZogSo, loga, v\ 
space. The thin line on which the galaxies are distributed in this space is the intersection of these 



S(s) = Sq exp 




Figure 1: Two by two correlations of the Sersic parameters logY,o,loga,i> for elliptical galaxies. Distribution 
of elliptical galaxies in the 3D Sersic parameter space (top right). The lines superimposed in each panel are 
the correlations predicted from the intersection of the Entropic Surface (constant specific entropy) and Potential 

energy-mass relation scaling relation. 



two two-manifolds (Fig. ||). Such relations are most probably a consequence of the formation 
and evolution processes undergone by these objects, since theory predicts U oc M 5 / 3 under the 
hypothesis that energy and mass are conserved (Marquez et al. 2001). 




Surface 

Figure 2: 3D representation of the specific entropy and potential energy-mass two-manifolds for elliptical galaxies, 

using the coordinates [log'Eo, loga, v\. 



Interestingly, numerical simulations of cold dark matter haloes in two different mass ranges 
lead to a similar scaling law between the potential energy and mass of the haloes. In the mass 
range 4 10 5 < M < 4 10 8 M & (unvirialized clusters), Jang Condell & Hernquist (2001) find a 
relation consistent with U oc M 5//3 , while in the mass range 10 12 < M < 10 15 Mq (virialized 
clusters), Lanzoni (2000) finds U oc M im±om . 

Since the X-ray gas in clusters of galaxies is weakly dissipative, it is likely to verify similar 
properties. The work presented here is meant to check this hypothesis. Our results are still 
preliminary and will be fully presented in a forthcoming paper (Demarco et al. in preparation). 



2 The data and data reduction 



We have selected a sample of 24 clusters with redshifts between 0.01 and 0.3 from the ROSAT 
PSPC archive with a high exposure time (and therefore signal to noise ratio), a roughly regular 
structure (obviously interacting clusters were discarded) and a proper light curve (no count rate 
peaks larger than 3 counts / seg in the 4 bands used) . The main properties of the selected clusters 
are given in Table |l|. 



Table 1: Physical parameters of the 24 clusters of our sample. 
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The data reduction was done using Snowden's now standard software (Snowden et al., 1994). 

3 The method 

3.1 Fitting the Sersic parameters 

In order to model the X-ray emission of the cluster, we have chosen a 3D density profile p(r), 
which is a generalized form of a Mellier-Mathez profile (Mellier &: Mathez 1987). It is completely 
determined by its three parameters and can be written under the form: 

p(r) = po (r/ay p exp[-(r/a) u ] (2) 
where the parameters p and v are related by the numerical aproximation (Marquez et al. 2001): 



p ~ 1.0 - 0.6097i/ + 0.05563i/ 2 



(3) 



This relation gives the best approximation to the Sersic law when equation (|2|) is projected. 
The Sersic profile defined by equation (|l]) corresponds to a surface mass density while equation 
(|2|) is the volume mass density. The condition that the mass obtained by integrating equation (|I]) 
must be equal to the mass obtained by integrating equation (0), gives us the following relation 
between the parameters So and po: 

1 r(2) 

where T(a) is the gamma function defined by T(a) = x a ~ 1 e~ x dx. 

In order to determine the correct set of values for the three Sersic parameters for each cluster, 
we fit the ROSAT images by a pixel-to-pixel method which: i) creates a three-dimensional model 
of the X-ray emission and ii) projects it by integration along the line of sight, taking into account 
the energy response and the point spread function of the detector. The result is a synthetical 
image which can be compared with the observation. 

To obtain this synthetical image, we have used a code that takes into account the generalized 
Mellier-Mathez density profile (equation @) and not only the free- free Bremsstrahlung emission, 
but also the free-bound and bound-bound X-ray emissions. The code computes the X-ray 
emissivity, e u , in every point of the space; this emissivity is then projected by integration along 
the line of sight to obtain a surface brightness: 

A*( s ) = / / e f{ s + z ) w(v) dis dz 

J Z= OO «/f m j n 

where s 2 = x 2 + y 2 , with x and y the perpendicular directions to the line of sight, and z the 
coordinate along the line of sight. This integral is computed taking into account the energy 
response of the satellite w(v). 

Finally, this projected image is convolved with the ROSAT point spread function (PSF), 
which varies as a function of position on the detector and energy. We have used a FWHM of 2 
pixels which corresponds to the central PSF, because the resulting spreading is more important 
in the center where the intensity gradient is stronger. 

The cluster redshift and the gas temperature are required as input parameters. The redshift 
for each cluster was taken from the SIMBAD data base (except for A2199 for which the redshift 
was obtained from Wu, Xue & Fang (1999), and the gas temperature was taken from Wu, Xue & 
Fang (1999), except for A2034 for which the temperature was taken from Ebeling et al. (1996), 
and for A2382 for which we have used a temperature of 5.0 keV . We assumed a standard CDM 
cosmology with a Hubble constant Hq = 50 km s _1 Mpc^ 1 , £1q = 1, and A = 0. 

To obtain an initial guess for the free parameters in equation (^) we fit a Sersic profile 
(equation (|l])) to the X-ray surface brightness of each cluster. The ellipticity and semi-major 
axis position angle of the X-ray gas were also given as starting values for the code. 

Once we have obtained the first synthetic image with the initial guess for the value of each 
free parameter, the code compares it to the actual ROSAT image. The values of the parameters 
are then changed and the comparison process is repeated iteratively, until it finds the 3D X-ray 
emission (density dependent) which best fits the surface brightness profile of the observation 
when projected. The fitting process is carried out with the MIGRAD methode in the MINUIT 
library of CERN (James 1994). 

The best fit parameters for each cluster are: the semi-major axis a (in kpc), the shape 
parameter v (dimensionless) and the central electronic number density, n e Q (cm -3 ). 

Since the formula used to calculate the specific entropy is valid only for spherical symmetry, 
and in order to take the ellipticity of the X-ray emission into account, we have also computed 



(4) 



a new scale parameter, a eq , denned by a eq = a \fbfa. In this way, we will calculate the spe- 
cific entropy of a spherically symmetric X -ray region with a scale parameter a e q. The best fit 
parameters are given in Table 1. 



3.2 Calculation of the various physical quantities 

By integrating equation (^) we find the gas mass as a function of radius: 

M gas (r) = f p{u)^u 2 du = ±^ 7 (iz£ |( !l n ( 5 ) 
Jo v v a 

where 7(0, z) is the incomplete gamma function defined by 7(0, z) = x a ~ 1 e~ x dx. 

In order to estimate the total gas mass, we only have to integrate the last equation until 
r = 00 (note that p(r) does not diverge). The result is: 

M gas = 2vra 2 S ir(-) (6) 
v v 

which depends only on the free parameters of the Sersic law. The potential energy of the gas is 
given by: 

TT — lvl gas A 2 3 v 2 \v \vJl f7\ 

U P ot = -5— = 4 vr a L — (7) 

Kg K g 

where we have used a gravitational constant G = 1 and where the gravitational radius R g is 
defined by R g = a R*. Here a is the scale parameter and R* is a dimensionless radius given by 
the numerical approximation: ln{R*) ~ °- 82032 ~°-9 2446 Mg) _|_ 0.84543 (Marquez et al. 2001). 

The second principle of thermodynamics states that a dynamical system in an equilibrium 
state has a maximum entropy. For a self-gravitating system in quasi-equilibrium, such as a 
galaxy cluster, there is no state with an entropy maximum. On the contrary, the entropy is 
constantly growing, but slowly, in a secular time scale, during which we can consider it as 
nearly constant. In spite of the X-ray emission, which is responsible in many cases for cooling 
flow processes which affect the equilibrium state of the cluster, we may therefore consider these 
objects as structures where dissipation processes are negligible compared to their gravitational 
energy, thus settling into a thermodynamical equilibrium. 

The specific entropy of the intra-cluster gas can be obtained from the distribution function 
in the phase space, /(if, if), of the gas particles. The calculation is made using the microscopic 
Boltzmann-Gibbs definition: 

S _ffln(f)d 3 xd 3 v 



M gas Jfd^xdH 

Estimating the distribution function / may be very difficult and some hypotheses are needed 
to simplify the problem. The first hypothesis is that our system has spherical symmetry, and the 
second one is that the velocity dispersion of the gas particles is isotropic. With these assumptions 
we are neglecting any possible rotation of the gas, and / can be obtained by Abel inversion from 
the density profile (Binney & Tremaine, 1987). Using equations (|2|) and (|j), the gas specific 
entropy (Lima Neto et al., 1999; Marquez et al., 2000) is then: 

3 9 

s = -ZnE +-lna +F(y) (9) 
where F is a function of the v parameter given by the numerical approximation: 

1 34 / 1 \ 129 2 

F(v) ~ -0.795 ln{y) - — + 3.85 - + ln(T[-\) - 0.822 

V \v J V 



4 Results 



. 1 Correlation of the Sersic parameters two by two 
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Figure 3: Correlation between the shape and scale parameters. 
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Figure 4: Correlation between the shape and intensity parameters. 
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Figure 5: Correlation between the scale and intensity parameters. 



We see in Figs. [3|, ||| and || that the three Sersic parameters are correlated two by two, as 
was previously observed for the optical surface brightness of elliptical galaxies. In particular, 
a clear correlation exists between the shape and scale parameters (Fig. |||), with a shape very 
similar to that found for elliptical galaxies. The other parameters are also correlated, although 
with more dispersion than for ellipticals. 



4-2 Constancy of the specific entropy 



Once the scale parameter a eq , the shape parameter v and the central density po are determined, 
we can compute the specific entropy s for the gas component of the cluster from equation (|9|). 
The gas mass can be computed from equation (0). We can then calculate the gas entropy 
S = sM gas . 

A plot of the gas entropy 5 as a function of the gas mass M gas is shown in Fig. ||; a linear 
relation exists between the entropy and total mass of gas, implying that the specific entropy s 
is constant. The slope gives the value of s = 34.5 ±1.3. A remarkable characteristic of Fig. |6| 
is the very low dispersion of the points. 
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Figure 6: Relation between the gas entropy and mass. 



4-3 Scaling relation between the cluster potential energy and the X-ray gas mass 

The potential energy U po t of the gas is displayed as a function of the gas mass in Fig. [?]. 

If we write this relation as log(U) — a log(M) = const, we find a slope a = 1.72 ±0.05, which 
is very close to the theoretical value of 5/3. 

This result confirms the hypothesis made for the formation of structures. The fact that 
elliptical galaxies and gas in clusters of galaxies verify the scaling law between potential energy 
and mass before mentioned is a strong argument in favour of the idea that self-gravitating 
structures are formed by processes where the conservation of energy and mass are verified. 

4-4 Second order relations 

Numerical simulations of elliptical galaxies formed in a hierarchical merging scheme show that 
the specific entropy slightly varies with mass, in a similar way as for observed galaxies (Marquez 
et al., 2000). This weak dependence has been interpreted as due to merging processes (Lima 
Neto et al, 1999). 

We have searched for a similar relation in our sample of clusters. The plot of the gas specific 
entropy s versus the gas mass M gas shown in Fig. || clearly reveals a linear relation between 
s and log(Mg as ). The gas specific entropy s is therefore a function of the gas mass. However, 
this is a second order relation compared to the more important relation s = constant previously 
found. The difference with elliptical galaxies is that the slope in Fig. || is steeper for clusters 
than for ellipticals. Thus, if we write s = sq + (3 log(M), we find (3 = 4.7 ± 0.3 for clusters and 
/9 ~ 1 for ellipticals. 
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Figure 7: Scaling law between gas the potential energy and mass. The solid line corresponds to the best fit to 
the data and the dashed line is the theoretical curve (slope 5/3 in logarithmic units). 
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Figure 8: Relation between the gas specific entropy and mass. The straight line corresponds to the best fit slope 

/3 = 4.7. 



5 Conclusions 

The similarity of the relations found for the optical light distribution in elliptical galaxies and for 
the X-ray gas emission in clusters seems to confirm the hypothesis that the physical conditions 
prevalent when these self-gravitating systems were formed are comparable. These self-gravitating 
and almost dissipationless systems are likely to have evolved in a comparable way. 

The more important dependence of the specific entropy with mass for clusters suggests that 
dissipating processes such as Bremsstrahlung emission (L oc M 2 ) and cooling flows may play 
an important role as generators of entropy. Merging processes between clusters may also be of 
some importance. 

The present work was carried out on a sample of 24 nearby clusters; the study of more distant 
objects could give us information on a possible dependence of specific entropy with redshift and 
on a possible evolution of the Sersic scale parameter a. Such an analysis obviously requires 
a higher spatial resolution and sensitivity than those of ROSAT, and should be possible with 
XMM or Chandra. 
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